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Abstract 



We study the general Higgs-weak boson coupling with CP-violation via the 
process e + e~ — > ffh. We categorize the signal channels by sub-processes Zh 
production and ZZ fusion and construct four CP asymmetries by exploiting 
polarized e + e~ beams. We find complementarity among the sub-processes and 
the asymmetries to probe the real and imaginary parts of the CP-violating 
form factor. Certain asymmetries with unpolarized beams can retain signif- 
icant sensitivity to the coupling. We conclude that at a linear collider with 
high luminosity, the CP-odd ZZh coupling may be sensitively probed via 
measurements of the asymmetries. 
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I. INTRODUCTION 



Searching for a Higgs boson has been a major motivation for many current and future 
collider experiments, since Higgs bosons encode the underlying physics of mass generation. 
In the minimal Standard Model (SM), there is only one CP-even scalar. In the two-Higgs 
doublet model or the supersymmetric extension of the SM, there are two CP-even states, one 
CP-odd state, plus a pair of charged Higgs bosons. The couplings of Higgs bosons to elec- 
troweak gauge bosons are particularly important since they faithfully represent the nature 
of the electroweak gauge symmetry breaking. Determining the detailed properties of the 
Higgs boson couplings will be of fundamental importance to fully construct the theoretical 
framework of the electroweak sector. 




FIG. 1. Vertex function for ZZh coupling. 

The most general interaction vertex for a generic Higgs boson (h) and a pair of Z bosons, 
Z^(ki) Z l/ (k 2 ) h, can be expressed by the following Lorentz structure 

r^fo, k 2 ) =i 2 -h\a Mlg» v + b {k»k u 2 - h ■ k 2 g^) + b e» vp °k lp k 2a \, (1) 

where v = {^/2G F )- 1 / 2 is the vacuum expectation value of the Higgs field, and the Z 
boson four-momenta are both incoming, as depicted in Fig. [TJ. The a and b terms are CP- 
even and the b term is CP-odd. Thus, the simultaneous existence of terms a (or b) and b 
would indicate CP violation for the ZZh coupling |TJ-|3]]. We note that in the SM at tree 
level, a = 1 and b = b = 0. In supersymmetric theories with CP-violating soft SUSY 
breaking terms f|, these CP- violating interactions may be generated by loop diagrams. 
More generally, the parameters can be momentum-dependent form factors and of complex 
values to account for the dispersive [Re(b)} and absorptive [2m (6)] effects from radiative 
corrections. Alternatively, in terms of an effective Lagrangian, the b term can be from 
gauge invariant dimension-6 operators || , and the b term can be constructed similarly with 
CP-odd operators involving the dual field tensors. Dimensional analysis implies that the 
parameters b and b may naturally be of the order of (v/A) 2 where A is the scale at which 
the physics responsible for the electroweak symmetry breaking sets in, presumably A < 4irv. 
The CP-odd coefficient b is of course very much model-dependent. 

Possible CP-violation effects via Higgs-gauge boson couplings have recently drawn a lot 
of attention in the literature. In Ref . [[J , CP-odd observables in decays h — > ZZ, W + W~ 
and tt were constructed. It was discussed extensively how to explore the Higgs properties 
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via the process e + e _ — > Z7t at future linear colliders. The polarized photon-photon 
collisions for 77 — * h and the electron-electron scattering process e~e~ — > e~e~h [[/J were 
also considered to extract the CP-violating couplings. There has also been considerable 
amount of work for investigation of CP- violating Higgs boson interactions with fermions at 
future e + e~ colliders ||. 

In this paper, we study the CP-violating coupling of ZZh at future e + e~ linear colliders. 
In Section II, we set out the general consideration, identifying the Zh production and ZZ 
fusion signals and exploring the generic CP-odd variables by exploiting the polarized beams. 
Given specific kinematics of the signal processes under investigation, we construct four CP 
asymmetries in Section III. We find important complementarity among the sub-processes 
and the asymmetries in probing different aspects of the CP-odd coupling, namely the real 
(dispersive) and imaginary (absorptive) parts of b. We also examine to what extent this 
coupling can be experimentally probed via measurements of the CP asymmetries, with 
and without beam polarization. We present some general discussions of our analyses and 
summarize our results in Section IV. 



II. GENERAL CONSIDERATION 

We concentrate on the scenario with a light Higgs boson below the VT^-pair threshold. 
The Higgs-weak boson coupling will be studied mainly via Higgs production, rather than its 
decay. We focus on the Higgs boson production associated with a fermion pair in the final 

state 

e-(pi)e + (p 2 )^f( qi ) f(q 2 )h(q 3 ). (2) 
The Higgs boson signal may be best identified by examining the recoil mass variable 

m rec = (Pi +P2-Q1- (I2) 2 = s + m 2 ff - 2y/s(E f + Ef), (3) 



where rriff is the // invariant mass and Ef {Ef) is the fermion (anti- fermion) energy in the 
c. m. frame. This recoil mass variable will yield a peak for the signal at the Higgs mass m h , 
independent of the Higgs decay. This provides a model-independent identification for the 
Higgs signal. For this purpose, we will accept only 

f = e~ , and u, d, s (4) 

to assure good energy determination for the final state leptons and light quark jets. When- 
ever appropriate, we adopt energy smearing according to a Gaussian distribution as 

AE 12% . , 

— — = -= © 1% for leptons (5) 
E ^E 

' ° ° © 2% for quarks. (6) 



E 

In realistic experimentation, the charged tracking information may also be used to help 
improve the momentum determination. 
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FIG. 2. Normalized mass distributions for e + e~ — > e + e~h at ^/s = 500 GeV with 
m h = 120 GeV. The dashed curve is for the recoil mass in Eq. (|3|), and the solid is for the 
invariant mass m ee . 



As an illustration, the recoil mass spectrum for an e + e~ final state is shown in Fig. |2| 
by the dashed curve. The width of the peak in m rec spectrum is determined by the energy 
resolution of the detector as simulated with Eq. (H). We have also required the final state 
fermions to be within the detector coverage, assumed to be 

| cos 6 f \ < cos 10° (7) 

with respect to the beam hole. 



A. Zh production versus ZZ fusion 

The signal channel Eq. (fj) can be approximately divided into two sub-processes 

e + e~ — > Zh (Zh production) , (8) 
Z*Z* h (ZZ fusion). (9) 

Eq. (H) yields light fermion states of all flavors from Z decay; while Eq. @ always has an 
e + e~ pair in the final state. These two sub-processes can be effectively distinguished by 
identifying the final state fermions. Even for the final state of e + e~, one can separate them 
by examining the mass spectrum m ee . This is illustrated in Fig. ^|for the e + e~ final state by 
the solid curve. The sharp peak at M z indicates the contribution from the decay Z — > e + e~, 
while the continuum spectrum at higher mass values is from ZZ fusion sub-process. In our 
analysis, we have included both contributions coherently. However, when necessary, we 
separate out the ZZ fusion contribution by requiring 
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FIG. 3. Total cross sections for e + e~ — > e + e~h in fb (a) versus for representative val- 
ues of m h , and (b) versus m h for representative values of \/s. The dashed curves are for 
e + e~ — ► Zh — > e + e~h only. No kinematical cuts are imposed. 

m ee > 100 GeV. (10) 

The Zh associated production is the leading channel for Higgs boson searches at e + e~ 
colliders. ZZ fusion, on the other hand, is often thought to be much smaller due to the 
small Zee vector coupling and low radiation rate of Z bosons off e beams. However, the 
rate of the fusion process increases with c. m. energy logarithmically like ln 2 (s/M|), and it 
is also more important for higher Higgs masses. The ZZ fusion process naturally leads to 
a pair of electrons in the final state, which is desirable when the charge information of the 
final state is needed. Moreover, due to the helicity conservation at high energies, the Zh 
production has only helicity combinations for the initial e + e~ of (H — ) and ( — h); while the 

ZZ fusion has ( ) and (++) in addition, where — (+) refers to the left (right) handed 

helicity. These additional helicity amplitudes may provide further information regarding the 
CP test, as we will see in the later analysis. 

Figure |3] presents the total cross sections for e + e~ — > e + e~h to demonstrate the compar- 
ison between the Zh and ZZ fusion processes. Figure |3](a) gives cross sections versus \/s for 
m h = 110 — 200 GeV, and (b) versus m h for s/s = 350—800 GeV. The solid curves are for the 
total SM rate including all contributions coherently, and the dashed curves are with a real Z 
decay for Zh -> e + e~h. We see that at y/s = 500 GeV and m h = 120 GeV, a(ZZ h) w 10 
fb> 2a(Zh -> e+e~/i, /i + /i~/i). At y/s = 800 GeV and m h = 120 GeV, the fusion cross 
section becomes about an order of magnitude higher than that of Zh —>■ e + e~h, [i + fi~h. 
Clearly, at a linear collider above the Zh threshold, the ZZ fusion process is increasingly 
more important in studying the Higgs properties 0. 
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B. CP property 



To unambiguously identify the effect of CP violation, one needs to construct a "CP-odd 
variable" , whose expectation value vanishes if CP is conserved . We begin our analysis 



by examining the CP-transformation property. First of all, we note that the initial state of 
Eq. (0) can be made a CP eigenstate, given the CP-transformation relation 

e~(ax,p) e + (>2, -p) e~(-a 2 ,p) e + (-cr 1 , -p), (11) 

where <jj is the fermion helicity. Now consider a helicity matrix element .M CT1CT2 (<fi, where 
G\ (a 2 ) denotes the helicity of the initial state electron (positron), which coincides with the 
longitudinal beam polarization; q\ (q 2 ) denotes the momentum of the final state fermion 
(anti- fermion). It is easy to show that under CP transformation, 

-M-+(gi, q 2 ) => M-+{-q 2 , -qi), (12) 
M—(qi,q 2 ) M++(-q 2 ,-qi), (13) 



and .M+-, M.++ transform similarly. If CP is conserved in the reaction, Relations ([12]) and 



(O) take equal signs. These relations precisely categorize two typical classes of CP test: 



CP eigen-process: Under CP, M. |_ (or M.+-) is invariant if CP is conserved. One can 

thus construct CP-odd kinematical variables to test the CP property of the theory. We can 
construct a "forward-backward" asymmetry 

A FB = a F -a B = ( l -^—d cos 9 - f° -^—d cos 9, (14) 
Jo a cost' J -1 a cost' 

with respect to a CP-odd angular variable 9 . This argument is applicable for unpolarized 
or transversely polarized beams as well. 

CP-conjugate process: M. and M. ++ are CP conjugate to each other. In this case, instead 

of a kinematical variable, the appropriate means to examine CP violation is to directly 
compare the rates of the conjugate processes. We can thus define another CP asymmetry 
in total cross section rates between the two conjugate processes of opposite helicities, called 
the "left-right" asymmetry 

A L r = o— - (T++. (15) 

The longitudinally polarized cross section for arbitrary beam polarizations can be calcu- 
lated by the helicity amplitudes 

da(P_P + ) = + P_)(l + P + )da ++ + (1 + P_)(l - P + )do + _ 

+ (1 - P_)(l + P + )Ax_+ + (1 - P_)(l - P+)tZ<7__], (16) 

where P_ (P+) is the electron (positron) longitudinal polarization, with P± = —1 (+1) for 
purely left (right) handed. Whenever appropriate in our later studies, we will assume the 
realistic beam polarization as (|P_|, \P+\) = (80%, 60%) Jll|. 



6 




FIG. 4. Normalized polar angle distributions for a |_ at \fs = 500 GeV with m h = 120 GeV 

for (a) e + e~ — > Zh with Z — * ff, and (b) e + e~ — > e + e~h via ZZ fusion. The solid curves are 
for the SM interaction (a = 1), the dashed for the CP-odd (Xm(b) = 1), and the dotted for CP 
violation with a = Xm(b) = 1. 100% longitudinal polarization of e^e^ has been used. 



III. CP-ODD VARIABLES AND THE CP-ODD COUPLING 



In this section, we construct CP-odd variables for the Higgs signal in Eq. (FJ) in order 
to study the CP- violating interactions in Eq. ([!]). Different CP asymmetries appear to be 
complementary in exploring different aspects of the CP-odd coupling b. 




A. Simple polar angles and Xm(b) 

It has been argued that the Zh process will test the spin-parity property of the 
coupling by simply measuring the polar angle distribution of the outgoing Z boson. The 
distribution can be written in the form 

scalar h 

d cos 9z \ 1 — | sin 2 6z pseudo-scalar h 

In fact, this simple polar angle may provide CP information as well. If we rewrite this angle 
in terms of a dot-product 

pi ■ q+ 

cos Z = | ^ | , (17) 

N \q+\ 

where q + = qi+cfe is the vector sum of the outgoing fermion momenta, it is easy to see that it 
is P-odd and C-even under transformation for the final state. One could thus expect to test 
CP property of the interactions by examining the polar angle distribution. The experimental 
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FIG. 5. Forward-backward asymmetries for a y versus 2m(b) at ^/s = 500 GeV with 

m h = 120 GeV for (a) Zh — > ffh: asymmetry in fb, (b) Zh —* ffh: percentage asymmetry, (c) 
ZZ fusion: asymmetry in fb, and (d) ZZ fusion: percentage asymmetry. The dashed curves are for 
100% longitudinal polarization e^e^>, the solid for a realistic polarization (e^,e^) = (80%, 60%), 
and the dotted for unpolarized beams. The error bars are statistical uncertainties obtained with a 
luminosity of 1000 fb -1 . 

study is made particularly simple since this variable does not require charge identification 
for the final state fermions. Because of this, one expects to increase the statistical accuracy 
by including some well-measured hadronic decay modes of Z, as we accept the light quark 
jets of Eq. ([|). However, after the azimuthal angle integration the dispersive part of the 
form factor proportional to 73e(5) vanishes and the surviving term is the absorptive part 
proportional to 2m(b). The angular distributions are shown in Fig. |](a) for ^Js = 500 GeV 
with m h = 120 GeV. The solid curve is for the SM interaction only (a = 1), the dashed 
curve is for the CP-odd only (Xm(b) = 1), the dotted is for CP violation with a = Im(b) = 1. 
We see from the dotted curve that there is indeed an asymmetry with respect to the forward 
(tt/2 < 9z < 0) and backward (n < 9z < vr/2) regions. We have assumed 100% longitudinal 
polarization of e^e^ for illustration here. 

The above calculation can in principle be carried through for the ZZ fusion process. 
However, due to the unique kinematics in this process, it appears that we can define an 
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FIG. 6. Polarized cross section asymmetry at y/s = 500 GeV with m h = 120 GeV versus Im(b) 
for (a) the asymmetry in fb, and (b) the percentage asymmetry. The dashed curves are for 100% 
beam polarization, the solid for a realistic polarization (e~,e + ) = (80%, 60%). The error bars are 
statistical uncertainties obtained with a luminosity of 1000 fb -1 . 



alternative polar angle 



cosfl_ = %^#^# , (18) 
|pi x q_\ ft x q 2 \ 



where g_ = q\ — q\, which yields a larger asymmetry and thus being more sensitive to 
the coefficient Tm(b). It is easy to verify that this variable is P-odd and C-even under 
transformation for the final state. Figure |](b) shows the angular distributions for e + e _ — > 
e + e~h via ZZ fusion with 100% longitudinal polarization of e^e\. The legend is the same 
as in (a). We see from the dotted curve that an asymmetry exists with respect to this angle. 

Replacing 9 by 6z in Eq. flU]), we can define a forward-backward asymmetry Aq B with 
respect to the angle 9z, and similarly another asymmetry Aq B with respect to the angle 
0_. These two asymmetries are calculated for cr_ + , and shown in Fig. |5| at yfs = 500 GeV 
with m h = 120 GeV versus Im(b). Figures |](a) and (b) are the asymmetry in fb and the 
percentage asymmetry respectively, with respect to 6z in Zh — > ffh. Similarly, Figs. |5](c) 
and (d) show the asymmetry and percent asymmetry for ZZ fusion with respect to 6L. 
The dashed curves are for 100% longitudinal polarization e^e^, the solid are for a realistic 
polarization (e£,eji) = (80%, 60%), and the dotted are for unpolarized beams. We see 
that the beam polarization here substantially enhances the asymmetries, and the realistic 
polarization maintains the asymmetries to a large extent. Some degree of asymmetry still 
exists even for unpolarized beams. The percentage asymmetry for the Zh process can be as 
large as 30% for Xm(b) ~ 0.2, and is typically of a few percent for ZZ fusion. 

We wish to address to what extent an asymmetry can be determined by experiments. 
For this purpose, we estimate the statistical uncertainties for the asymmetry measurements. 
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FIG. 7. Cross section asymmetries in fb versus y/s with m h = 120 GeV for (a) for- 
ward-backward asymmetry with respect to 9z for 2m(b) = 0.1, (b) forward-backward asymmetry 

with respect to 9- for for Xm(b) = 0.5, (c) LR asymmetry between a and a — h+ for Zm{b) = 0.1. 

The dashed curves are for 100% longitudinal polarization, and the solid for a realistic polarization 
(e~,e + ) = (80%, 60%). The error bars are statistical uncertainties obtained with a luminosity of 
1000 fb -1 . 



We determine the Gaussian statistical error by \J Np + Nb where Np (Nb) is the number 
of forward (backward) events. The statistical significance for the asymmetry measurement 
is obtained by 



\N F -N B \ 



(19) 



The error bars in the plots are calculted with an assumed integrated luminosity of 1000 fb -1 . 
Due to the larger asymmetry as well as a larger cross section for Zh — > ffh, the Zh 
production would provide a much better determination of Im(b). 

As we discussed earlier, the ZZ fusion process can provide another type of asymmetry 

between CP conjugate processes, in particular between a and a ++ as denned in Eq. (|T5D, 

which is absent in Zh production. This is presented in Fig. || for Alr, at y^i = 500 GeV 
with m h = 120 GeV versus Tm(b). Figure ^|(a) is the asymmetry in fb. The legend is the 
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FIG. 8. Normalized angular distributions for a |_ at = 500 GeV with m h = 120 GeV for 

(a) e + e~ — > Zh with Z — ► e~e + ,fi~^ + , and (b) e + e~ — ► e + e~h via ZZ fusion. The solid curves 
are for the SM interaction (a = 1), the dashed for the CP-odd (73e(6) = 1), and the dotted for CP 
violation with a = 'Re(b) = 1. 100% longitudinal polarization of e^e^ has been used. 

same as in Fig. |5|. The error bars are for a total integrated luminosity of 1000 fb _1 (500 fb _1 

each for a and cr ++ ). The percentage asymmetry in Fig. ||(b) can be at a 10% level for 

Xm(b) ~ 0.2. It is interesting to note that the solid curves yield a non-zero value for 6 = 0. 
This is due to the intrinsic LR asymmetry of the Z coupling to electrons. This shift appears 

when |P_| 7^ \P + \ and is proportional to a |_ — o~+_. It can be well predicted in the SM for 

a given beam polarization. 

Cross section asymmetries versus y/s are shown in Fig. [7] in units of fb with m h = 
120 GeV (a) forward-backward asymmetry for <t_ + in Zh production with respect to 6z for 
Tm{b) =0.1, (b) forward-backward asymmetry for cr_ + in ZZ fusion with respect to 6L for 

for Tm(b) = 0.5, (c) LR asymmetry between a and a ++ for Xm(b) = 0.1. The dashed 

curves are for 100% longitudinal polarization, and the solid are for a realistic polarization 
(e~, e + ) = (80%, 60%). The error bars are for the statistical uncertainty with a luminosity of 
1000 fb -1 . We see again the possibly good accuracy for determining the asymmetry by the 
Zh process. Furthermore, these two processes are complementary: at lower energies near 
threshold the Zh production is far more important; while at higher energies the ZZ fusion 
becomes increasingly significant, as has been seen in Fig. [| In Fig. 0(c), the reason that the 
realistic asymmetry (solid) is even bigger than the ideal case (dashed) is due to the non-zero 
contribution from the CP-conserving LR asymmetry of the Z coupling as discussed in the 
last paragraph. 
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FIG. 9. Forward- backward asymmetry for a y versus 72e(6) at = 500 GeV with 

m h = 120 GeV for (a) Zh — > ffh: asymmetry in fb, (b) Zh — * ffh: percentage asymmetry, (c) 
ZZ fusion: asymmetry in fb, and (d) ZZ fusion: percentage asymmetry. The dashed curves are for 
100% longitudinal polarization e^e^, the solid for a realistic polarization (e^,e^) = (80%, 60%), 
and the dotted for unpolarized beams. The error bars are statistical uncertainties obtained with a 
luminosity of 1000 fb -1 . 



B. The lepton momentum orientation and ~Re(b) 

We showed in the last section that the simple polar angles can probe CP violation for a 
Higgs-gauge boson coupling, but only for the absorptive part of the form factor Tm(b). In 
order to be sensitive to the dispersive part Ve(b), one needs to construct more sophisticated 
variables, involving the azimuthal angle information for the final state fermions. We find 
that a simple variable to serve this purpose |T2| can be defined as 



cos 9f 



Pi ■ %- x gg) 
l<?i x Qz\ 



(20) 



where q\ x q 2 defines the orientation of the plane for the final state fermion pair. This 
variable is P-even and C-odd under final state transformation. However, we would need 
to unambiguously identify the fermion from the anti-fermion, and to accurately determine 
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FIG. 10. Forward-backward cross section asymmetries for a |_ with respect to 0£ in fb ver- 
sus \fs with m h = 120 GeV and ~Re(b) = 0.1 for (a) Zh production and (b) ZZ fusion. The 
dashed curves are for 100% longitudinal polarization, and the solid for a realistic polarization 
(e~,e + ) = (80%, 60%). The error bars are statistical uncertainties obtained with a luminosity of 
1000 fb -1 . 



their momenta. This is naturally achievable for the ZZ fusion process, while we will have 
to limit ourself to / = e~,fi~ for the Zh — > ffh process. Explicit calculations show that 
this variable is only sensitive to 7?e(b) and insensitive to Tm(b). 

We evaluate the angular distribution for cos#£ at \fs = 500 GeV with m h = 120 GeV. 
Shown in Fig. || are the normalized distributions for (a) e + e~ — > Zh with Z — > e~e + ,fi~[i + , 
and (b) e + e~ — > e + e~h via ZZ fusion. The solid curves are for the SM interaction (a = 1), 
the dashed curves are for the CP-odd (73e(6) = 1), the dotted are for CP violation with 
a = TZe{b) = 1. 100% longitudinal polarization of e^e^ has been used as for cr_ + . The CP 
asymmetries are manifest as seen from the dotted curves. We define a CP asymmetry Ag B in 
the same way as in Eq. (0). The asymmetries for these two processes are calculated for a_ + , 
and shown in Fig. |9| at y/s = 500 GeV with m h = 120 GeV versus Ve(b). The parameters 
and legend are the same as in Fig. ||. We see that the percentage asymmetry for the Zh 
process is about 10% percent and for ZZ fusion it can be as large as 30% for Ve{b) ~ 0.2. 
The error bars in the plots are estimated with an integrated luminosity of 1000 fb -1 . Due 
to the large asymmetries both Zh production and ZZ fusion processes could provide a good 
probe to the coupling We(b). A particularly important result as indicated in Figs. H(c) and 
(d) is that the asymmetry for the ZZ fusion is rather insensitive to the beam polarization. 

Forward-backward cross section asymmetries for a_ + with respect to Qg are shown versus 
\J~s in Fig. [10] with m h = 120 GeV and Ve(b) = 0.1. Figure |T0"|(a) is the asymmetry for 
Zh production, and (b) is for ZZ fusion. We see again good sensitivity for measuring the 
asymmetry especially by the ZZ fusion process and at higher energies, which appears to 
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y/s (GeV) 500 500 800 800 
C (fb' 1 ) 500 1000 500 1000 

A^(Zh) [-+} 0.0028 0.0022 0.0043 0.0032 
lm(b) AfBjZh) [unpol.] 0.019 0.013 0.025 0.019 

a( b (ZZ) [-+] 021 016 019 0.13 
A L r(ZZ) 0.071 0.045 0.065 0.041 

Ajj'fjZh) [-+] 0.023 0.018 0.019 0.014 

7&(6) A%f(ZZ) [-+] 0.021 0.017 0.014 0.009 

Af B (ZZ) [unpol.] 0.024 0.018 0.016 0.010 



TABLE I. 95% C.L. limits on b from the CP asymmetries defined in the text at 
= 500, 800 GeV with m h = 120 GeV, for two representative luminosities C = 500, 1000 fb _1 . 
Realistic polarizations of (80%, 60%) are used unless specified as "unpolarized" . 

have very little dependence on the beam polarization. 

To further assess the linear collider sensitivity to b, we compare all the CP asymmetries 
and present in Table | the 95% Confidence Level (2a) sensitivity limits with m h = 120 GeV 
for two collider energies y/s = 500, 800 GeV and two choices of integrated luminosity 
C = 500, 1000 fb -1 . Realistic polarizations of (80%, 60%) are used unless specified for no 
beam polarization by "unpolarized" . We see that at a 500 GeV linear collider with a total 
luminosity of 1000 fb _1 , the CP-odd coupling form factor may be sensitively probed to a 
value of about Xm{b) ~ 0.0022 and 7te(b) ~ 0.017 at a 95% C.L. The coupling may even be 
probed without a beam polarization to a level of about Tm(b) « 0.013 and 7Ze(b) ~ 0.018. 
The beam polarization improves the sensitivity to Tm(b) by about a factor of 5 — 6 via 
Aj*(Zh), but does little to 7&(6) through A F e B (ZZ). At ^ = 800 GeV, the sensitivity in 
Zh process is slightly degraded. On the other hand, the sensitivity in ZZ fusion is enhanced 
by about a factor of two due to the larger cross section and larger asymmetry at higher 
energies. 



IV. DISCUSSIONS AND CONCLUSIONS 

Before summarizing our results, a few remarks are in order. First, in previous studies of 
the Zh process a common variable is defined as 

a %- x g+) • %- x gg) 

cosy + = -— — -, (21) 

|pi x q_\ | ft x q 2 \ 

where q + = ft + q\ = pz- This variable seems quite suitable for the Zh production since 
it is the azimuthal angle formed between the Zh production plane and the decay plane of 
Z — > // if the Z momentum is chosen to define the rotational axis. However, this variable is 
P-even and C-even under final state transformation and thus cannot provide an unambiguous 
measure for CP violation alone. One would have to analyze other angular distributions to 
extract the CP property of the interaction. 
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As a second remark, one may consider our analysis for the ZZ fusion similar to that 
in e~e~ collisions @, since the only tree-level Higgs boson production at e~e~ colliders is 
via the ZZ fusion mechanism |L3|]. However, an e~e~ initial state cannot be made a CP 



eigenstate as evident from the discussion of Eq. ( jlT|) . The explicit CP asymmetry in e~e~ 
collisions would have to be constructed in comparison with the conjugate e + e + reactions. 

Finally, although the ZZh coupling under current investigation is arguably the most 
important interaction in the light of electroweak symmetry breaking, other interaction ver- 
tices such as Z-fh and 77/i may be equally possible to contain CP violation induced by loop 
effects. Although the CP asymmetries constructed in this paper should be generically appli- 
cable to the other cases as well, we choose not to include those coupling in our analyses for 
the sake of simplicity. However, in terms of our ZZ fusion study, since the photon-induced 
processes 77 — > h, 7Z — > h would mainly give collinear electrons along the beams, our kine- 
matical requirement to tag e + e~ final state at a large angle will effectively single out the 
ZZh contribution. 

To summarize our analyses of possible CP violation for the interaction vertex ZZh, we 
classified the signal channel into two categories as Zh production with Z —*■ ff and ZZ 
fusion. We proposed four simple CP-asymmetric variables 



FB 
6z 
FB 



Ai B 



for Zh production, 
for ZZ fusion, 
for ZZ fusion only, 
for both Zh, ZZ. 



We found them complementary in probing the CP-odd coupling form factor b. The first 
three are sensitive to Tm(b), while the last one sensitive to Ve(b). A$ B yields the largest 
asymmetry for Xm{b) (see Fig. §), while Ag e B is the largest for lZe(b) (see Fig. |[), both 
reaching about 30% for \b\ ~ 0.2. The ultimate sensitivity to b depends on both the size 
of asymmetry and the signal production rate. As illustrated in Table |, at a 500 GeV 
linear collider with a total luminosity of 1000 fb , the CP-odd coupling may be sensitively 
probed to a value of about Tm{b) ~ 0.0022 and Ve(b) ~ 0.017 at a 95% C.L. with the beam 
polarization (80%, 60%). The coupling may even be probed without beam polarization 
to a level of about Tm(b) ~ 0.013 and 7?e(b) ~ 0.018. At a higher energy collider with 
\fs = 800 GeV, the sensitivity in Zh process is slightly degraded but that in ZZ fusion is 
enhanced by about a factor of two. 

Acknowledgments: We thank R. Sobey for his early participation of this project. We 
would also like to thank K. Hagiwara, W.-Y. Keung, G. Valencia and P. Zerwas for helpful 
discussions. This work was supported in part by a DOE grant No. DE-FG02-95ER40896 
and in part by the Wisconsin Alumni Research Foundation. 



15 



REFERENCES 



[1] D. Chang, W.-Y. Keung and I. Phillips, Phys. Rev. D48, 3225 (1993); A. Soni and 
R.M. Xu, Phys. Rev. D48, 5259 (1993); A. Skjold and P. Osland, Phys. Lett. B329, 
305 (1994). 

[2] V. Barger, K. Cheung, A. Djouadi, B.A. Kniehl, and P.M. Zerwas, Phys. Rev. D49, 79 

(1994); M. Kramer, J. Kuhn, M.L. Stong and P.M. Zerwas, Z. Phys. C64, 21 (1994). 
[3] K. Hagiwara and M. Stong, Z. Phys. C62, 99 (1994); J.P Ma and B.H.J. McKellar, 

Phys. Rev. D52, 22 (1995); A. Skjold and P. Osland, Nucl. Phys. B453, 3 (1995); W. 

Kilian, M. Kramer and P.M. Zerwas, Phys. Lett. B381, 243 (1996); K. Hagiwara, S. 

Ishihara, J. Kamoshita and B.A. Kniehl, Eur. Phys. J. C14, 457 (2000). 
[4] A. Pilaftsis, Phys. Lett. B435, 88 (1998); Phys. Rev. D58, 096010 (1998); A. Pilaftsis 

and C.E.M. Wagner, Nucl. Phys. B553, 3 (1999); S.Y. Choi, M. Drees and J.-S. Lee, 

Phys. Lett. B481, 57 (2000). 
[5] K. Hagiwara, S. Ishhara, R. Szalapski, and D. Zeppenfeld, Phys. Rev. D48, 2182 (1993). 
[6] B. Grzadkowski and J.F. Gunion, Phys. Lett. B294, 361 (1992); S.Y. Choi, K. Hagiwara 

and M.S. Baek, Phys. Rev. D54, 6703 (1996); G.J. Gounaris and G.P Tsirigoti, Phys. 

Rev. D56, 3030 (1997), Erratum-ibid. D58, 059901 (1998). 
[7] C.A. Boe, O.M. Ogreid, P. Osland and J.-Z. Zhang, Eur. Phys. J. C9, 413 (1999). 
[8] X.-G. He, J.P. Ma and B. McKellar, Phys. Rev. D49, 4548 (1994); B. Grzadkowski and 

J.F. Gunion, Phys. Lett. B350, 218 (1995); J.F. Gunion, B. Grzadkowski and X.-G. He, 

Phys. Rev. Lett. 77, 5172 (1996); S. Bar-Shalom, D. Atwood, G. Eilam, R.R. Mendel 

and A. Soni, Phys. Rev. D53, 1162 (1996); S. Bar-Shalom, D. Atwood and A. Soni, 

Phys. Lett. B419, 340 (1998); K.S. Babu, C. Kolda, J. March-Russell and F. Wilczek, 

Phys. Rev. D59, 016004 (1999); B. Grzadkowski, J.F. Gunion and J. Kalinowski, Phys. 

Rev. D60, 075011 (1999); Phys. Lett. B480, 287 (2000); E. Asakawa, S.Y. Choi, K. 

Hagiwara and J.S. Lee |hep-ph/0005313| . 
[9] J.F. Gunion, T. Han and R. Sobey, Phys. Lett. B429, 79 (1998). 
[10] G. Valencia, TASI 1994: 0235-270 (QCD161 : T45 : 1994), |hep-ph/941144l . 
[11] P. Zerwas, plenary talk presented at The 5th International Linear Collider Workshop, 

Fermilab, Oct. 2000. 

[12] The same variable has been used for different processes, see, e. g., A. A. Likhoded, G. 

Valencia and O.P Yushchenko, Phys. Rev. D57, 2974 (1998). 
[13] V. Barger, J. Beacom, K. Cheung and T. Han, Phys. Rev. D50, 6704 (1994). T. Han, 

Int. J. Mod. Phys. All, 1541 (1996). 



16 



